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ABSTRACT 

With 16-month Kepler data, 14 long-period (40 d — 265 d) eclipsing binaries 
on highly eccentric orbits (minimum e between 0.5 and 0.85) are recognized from 
their closely separated primary and secondary eclipses (At i n = 3 d — 10 d). These 
systems confirm the existence of a previously hinted binary population situated 
near a constant angular momentum track at P(l — e 2 ) 3 / 2 ~ 15 d, close to the tidal 
circularization period P C i rc . They may be presently migrating due to tidal dissi- 
pation and form a steady-state stream (~ 1% of stars) feeding the close-binary 
population (few-% of stars). If so, future Kepler data releases will reveal a grow- 
^ ing number (dozens) of systems at longer periods, following dN/ d lg P oc P 1 / 3 

with increasing eccentricities reaching e — > 0.98 for P — > 1000 d. Radial-velocity 
follow up of long-period eclipsing binaries with no secondary eclipses could offer 
a significantly larger sample. Orders of magnitude more (hundreds) may re- 
veal their presence from periodic "eccentricity pulses", such as tidal ellipsoidal 
variations, near pericenter passages. Several new few-day-long eccentricity-pulse 
candidates with long period (P=25 d — 80 d) are reported. 

Subject headings: 



1. Introduction 

How do close binaries (P < 10 d) form? An intriguing clue comes from their surprisingly 
high fraction having a tertiary companion, reaching as much as ~ 96% for systems with 
P < 3 d (Tokovinin et al. 2006). This lends support to the "Kozai Cycles plus Tidal Friction" 
(KCTF) formation channel advocated by several authors (Kiseleva et al. 1998; Eggleton & 
Kiseleva-Eggleton 2001; Fabrycky & Tremaine 2007). In KCTF, a tertiary companion on a 
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highly inclined outer orbit can excite the inner binary orbit, initially having a long period, 
to very high eccentricity during Kozai-Lidov cycles 1 , which allows efficient tidal dissipation 
to take place during close pericenter passages, eventually shrinking the orbital period to a 
few days. 

One clear implication of KCTF is that there should exist a population of highly-eccentric 
binaries in the process of migrating to their final close-in orbits. Such a prediction was 
studied in the context of a similar formation process of "hot Jupiters" (close-in Jovian-mass 
companions) by Socrates et al. (2011). In this Letter, we explore observational implications 
of such a population which is possibly responsible for close-binary formation. 

2. Steady-State Distribution 

In this section, we briefly review the features of a steady-state distribution of migrating 
binaries due to tidal dissipation in analogy to migrating Jupiters discussed in Socrates et 
al. (2011). Stars have been forming at a roughly constant rate over the age of the Galaxy, 
therefore one would expect a steady-state source of new distant binaries that constantly 
migrate to lower periods, "feeding" the close-binary population. Tidal dissipation drains 
orbital energy but hardly changes the angular momentum of an eccentric binary. While 
losing energy due to tidal dissipation, eccentric binaries migrate on tracks of constant orbital 
angular momentum corresponding to constant P(l — e 2 ) 3 / 2 = Pf. Pf equals the period of the 
final circularized close binary (once circular, tidal dissipation ceases to operate in the tidally 
locked system). The rate of energy dissipation is a strong function of binary separation, and 
for an eccentric orbit, most of the energy loss occurs during the pericenter passage. At high 
e, due to the constant angular momentum, the pericenter passages have approximately the 
same parabolic shape, and the energy lost each orbit, AE, is approximately the same. This 
implies an average energy loss rate E oc 1/P, which results in a steady state distribution 
dTV/dlgP oc E/E oc P 1/3 , where E oc P~ 2/3 . This distribution at high e is relatively 
independent of tidal dissipation physics. At moderate and low e, the distribution is sensitive 
to the poorly understood tidal dissipation processes. The distribution obtained for the 
commonly assumed equilibrium tide model (e.g, Hut 1981) is derived in Socrates et al. 
(2011). An additional complication is that angular momentum oscillates while the Kozai 
mechanism operates. As shown in Socrates et al. (2011), this does not significantly affect 



1 Recent studies on Kozai-Lidov cycles show that the parameter space that allows excitation to very 
high eccentricity is substantially wider than previously thought due to the commonly neglected octupolc 
contribution to the tertiary potential (Ford et al. 2000; Naoz et al. 2011; Katz et al. 2011; Lithwick & Naoz 
2011). 
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the resulting distribution as long as we consider narrow bins in angular momentum (or Pf) 
in which most of the dissipation occurs. In this Letter, we focus on the long-period, highly 
eccentric members in this population, which are simple to analyze theoretically and, as shown 
below, are easily distinguishable observationally. 

3. Past Observational Hints 

Making meaningful statistical studies of the binary distribution requires an un-biased 
sample. The commonly adopted study on multiplicity of solar-type primaries is from a 
complete, volume-limited (within 22 pc and spectral type F7-G9) sample of ~ 164 systems 
by Duquennoy & Mayor (1991). They identify three period regimes. Binaries with periods 
shorter than about 10 d have eccentricities close to zero and are thought to be efficiently 
circularized by the strong tidal dissipation at small separations. Those with periods longer 
than ~ 1000 d are found to be consistent with an Ambartsumian distribution, dN/de = 2e, 
which is expected if energy conserving dynamical perturbations operate. In the intermediate 
period range 10 d < P < 1000 d, the distribution in e is described as bell-shaped with an 
average of e ~ 0.3, which has since been widely accepted as reflecting the initial distribution 
at binary formation. They note that the largest eccentricities in the 10 d < P < 1000 d regime 
seem to be curiously populated by triple systems. They further note that just beyond the 
selection criterion of their complete sample, there exists some clear outliers to the eccentricity 
distribution, and interestingly two of them, HD 137763 (e = 0.975, P = 890 d) and HD 
114260 (e = 0.56, P = 20.5 d) are found to have additional companions. 

Raghavan et al. (2010) provide a new volume-limited sample selected using Hipparcos 
parallax (within 25 pc and stellar type F6-K3) including 454 solar-type stars, a factor of 
~ 3 larger in size than that of Duquennoy & Mayor (1991). Figure 1 shows the period 
and eccentricity distribution of binaries in this sample. Systems with known companions 
(multiples) are shown in filled circles while those without are plotted in crosses. As pointed 
out by the authors, the upper envelope of eccentricity for systems with 10 d < P < 1000 d 
appears to mostly consist of multiples, indicative of Kozai oscillations. The five "eccentric 
envelope" systems, each with eccentricity larger than that of any binary with shorter period, 
have the following parameters: (P, e, P te rtiary) = (890 d, 0.975, 2.4x 10 4 yr), (88 d, 0.81, 700 yr), 
(48 d, 0.61, 100 yr), (24 d, 0.41, 2000 yr), (21 d, 0.25, 260 yr) 2 . 

These envelope multiples, which possibly consist of ~ 1% of all main-sequence stellar 
systems, are very interesting in the context of KCTF. Figure 1 shows constant angular 



2 Orbital parameters taken from the 20 Apr 2010 updated version of Tokovinin (1997) 
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momentum tracks, as discussed in section 2, of Pf = P(l — e 2 ) 3 / 2 = 10 day and 30 day, 
which enclose the envelope. This raises the exciting possibility that these systems may 
represent a population of highly eccentric binaries in steady-state KCTF migration. If true, 
the fact that the Pf values of these systems are close to the circularization threshold of 
-Pcirc ~ 10 d is not surprising. Systems with Pf <C P C i rc , circularize quickly, while those with 
Pf ^> Pcirc, do not significantly migrate. Systems at Pf ~ P circ migrate on a time scale 
comparable to the age of the Galaxy, implying that a significant fraction of them should be 
currently migrating. 

Curiously, the tertiary companions of the above five systems all share Very long, Ptertiary ^ 
100 — 10 4 yr, orbital periods. More detailed examination of their orbital parameters indi- 
cate that due to the large perturber separations, the amplitude of the Kozai oscillations (if 
present) is suppressed by relativistic precession (see the discussions on "quenched" Kozai os- 
cillation in Socrates et al. 2011 and Figure 1 of Fabrycky & Tremaine 2007 for an example), 
which would naturally explain why they form an upper eccentric envelope. 

Unfortunately, the small-number (~ 5) statistics makes it difficult to robustly estab- 
lish the existence of this interesting eccentric envelope population and to investigate its 
properties. 

4. Kepler Probes 

Kepler mission has great potential to provide a complete sample that is orders-of- 
magnitude larger than Raghavan et al. (2010) to probe the above-mentioned population 
between 20 d and 1000 d that consists of possibly ~ 1% of main-sequence stellar systems. 
Kepler is a high-cadence, high-precision mission that aims to find habitable Earth-size plan- 
ets. It has been monitoring ~ 150, 000 FGK stars continuously with high cadence and high 
photometric precision since 2009, and is designed to operate for 3.5 yr with the possibility of 
an extension of several more years. The latest data public release available on 10 January, 
2012 includes 16-month long light curves. There are potentially ~ 1500 systems in total 
belonging to the eccentric envelope population among Kepler targets. 

4.1. Eclipsing Binaries 

A straightforward method to detect binary stars is to look for eclipses i.e., eclipsing bi- 
naries (EBs). For an isotropic distribution of orbital orientations, the chance of seeing an EB 
with at least one of the stars being eclipsed is Peb = Ptot[a(l — e 2 )]~ 1 (l + 2e/n) if the sum of 
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the radii of the two stars i? to t is much smaller than their separation. This geometrical prob- 
ability is roughly constant for a given orbital angular momentum J oc [a(l — e 2 )] 1 / 2 , and for 
a track with Pf = 15 day and high e, is approximately Veb = 0.06(i? to t/P©)(M;ot/^©) _173 ) 
where M tot is the total binary mass. We therefore expect to find of order ~ 100 EBs within 
the eccentric envelope population with 20 d < P < 1000 d. 

While the period is a direct observable of an EB light curve, the eccentricity is not. 
One way to accurately determine e is to make radial velocity observations. A second option 
is to carefully model the eclipse light curve using detailed characterization of the stellar 
parameters, allowing for a reasonable estimate of the eccentricity. In this letter we follow 
a third path by focusing on a subsample of EBs with both primary and secondary eclipses 
observed, upon which one can place strict lower limits on the eccentricity by measuring 
the timing separation between the two eclipses. We show that a robust selection of highly 
eccentric binaries can be made by searching for long period eclipses with short primary- 
secondary separations. 

The probability of observing both eclipses is given by Veb,2 = (tt — 2e)/(7r + 2e)"P EB , 
which amounts to about ~ 20% of all eclipsing systems at high e. We thus expect to find 
about two dozen systems from the eccentric envelope population . The timing difference 
between the eclipses, Atiji = tu — h is given by 



P 

Ati, n (P,e,u) = - 

7T 



ecoso; \ n K ecosw 



arecos ! . ! — VY— e 2 - 



V 7 ! - e 2 sin 2 u J ' ' 1 - e 2 sin w 



(1) 



where u> is the argument of pericenter. At high e the eclipses are most likely to occur near 
pericenter and the separation between the two eclipses is of order the pericenter passage time 
~ P(l — e) 3 / 2 /2. The separation reaches minimum at u = and is given by 

p 2P f 

Ati >n>min (P, e) = — (arccose - ey/l - e 2 ) > — J -. (2) 

and is approximated by P(l — e) 3 / 2 /2 to an accuracy of better than 20%. 

We search for EBs with periods longer than P cut = 40 d with secondary eclipses separated 
by less than Ati n, C ut = 10 d from the primary eclipses. Any EBs that survive the latter cut 
must have an eccentricity larger than e cut , which satisfies Ati i n imin (e cu t) = A£i ; n ;CUt . The 
resulting minimal eccentricity at Ati 5 n iCUt = 10 d as a function of P is shown in a curved 
dashed line in Figure 1. The period cut P cut = 40 d is also shown as a dashed line. Any 
observed EB that satisfies these observational cuts must reside in the (P, e) domain defined 
by these two lines, which allows for a robust selection of the eccentric envelope population . 
Moreover, > 80% of systems with Pf = 15 d and P > P cut = 40 d that show two eclipses, 
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satisfy At^n < Ati : ii jCU t = 10 d implying that this cut does not significantly affect systems 
in the eccentric envelope population . 

We search for EBs in the Janurary 2012 release of Kepler data, with all long-cadence 
(29.4min), Pre-Search data Conditioning (PDC) light curves from quarters Q1-Q6 (spanning 
16 months in total). First we identify eclipse candidates using a simple criterion that at 
least 3 consecutive points are above 3a deviations from the local median (determined with 
a window of 200 points). Then candidate light curves are selected by requiring pairs of 
candidate eclipses to repeat with periods longer than 40 d and that at least one eclipse be 
deeper than 1%. By visually inspecting each of the selected light curves, we identified 13 
binary systems with double eclipses having timing separations smaller than 10 d (and periods 
greater than 40 d). Six out of these are listed in the Kepler EB catalog (Slawson et al. 2011), 
which is based on the Kepler Q0-Q2 data (with a total time span of 125 day). Among the 
binaries in the catalog, we found an additional system that survives our period and eclipse 
separation cuts, but for which both eclipses have depths smaller than 1%. Therefore our 
selection appears to be efficient. The parameters of the 14 systems are listed in Table 1 and 
plotted as red arrows in Figure 1 showing their periods and minimal eccentricities derived 
from Equation (2), based on their measured At\ u. The light curves of the three systems with 
the longest periods are shown in Figure 2 for illustration. Clearly the search is incomplete at 
long periods, P > 160 d, at which fewer than 3 primary eclipses occur during the observation 
time span. Our detections establish that at least ~ 1% of stellar systems are binaries (with 
additional possible companions) in the selected region in the (P, e) plane. 



4.2. Eccentric Ellipsoidal Pulses 

Kepler's exquisite photometric precision opens up a venue of binary star detections 
with ellipsoidal, reflection/irradiation, and relativistic beaming modulations (e.g. Loeb & 
Gaudi 2003; Faigler et al. 2011). For very eccentric binaries, these modulations show up as 
short "pulses" that reach maximum amplitude near pericenter and last about the pericenter 
passage time. Welsh et al. (2011) reported the first Kepler detection of such "eccentric 
pulses" in KOI-54, with P = 41.8 d and e = 0.83. Thompson (2011) reported discoveries of 
17 systems with P < 20 d, dubbed as "heartbeat stars". 

In the case that ellipsoidal variations dominate (for small pericenter separation r p ), 
orbital parameters including eccentricity, argument of pericenter and the inclination can be 
robustly extracted from the shape of the pulses (Kumar et al. 1995). In contrast to EBs, the 
eccentric ellipsoidal variations can be virtually seen at all inclinations with similar amplitude, 
thus in principle allowing the detection of orders of magnitude more such systems than from 
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EBs. The amplitude of the ellipsoidal pulse due to a Sun-like primary in a long-period, highly 
eccentric system with mass ratio q and Pf ~ 15 d is roughly ~ q{Ro/r p ) 3 ~ 5 x 10 _4 g/(l +g) 
with a width of ~ 3d. Jenkins et al. (2010) demonstrates that at Kepler magnitude m^ cp < 
14 (half of all Kepler targets), the majority of the Kepler dwarfs reach better than ~ 10~ 4 
precision over 6.5 hr interval, including photometric noise and stellar variability. Therefore 
the prospect of detecting a large number of Sun-like binaries exhibiting eccentric ellipsoidal 
variations is promising. 

One simple way to find eccentric pulses is to examine the EB light curves (with either 
one or both of the primary and secondary eclipses) and inspect the region in the vicinity 
of the eclipses where such pulses are most likely to occur. The known periodicity helps 
distinguishing the pulses from other sources of variabilities. Figure 3 shows the folded light 
curve of the EB KIC 6864859 with detection of pulses with amplitude of 0.00025. The 
pulse shape, shown in the bottom panel, is well described by the linear-perturbation model 
of ellipsoidal variations 5F oc [1 — 3 sin 2 i cos 2 (/ + co)]/r 3 (Kumar et al. 1995), where / is 
the true anomaly and r is the instantaneous binary separation, for parameters (e = 0.63, 
bj = 2° and % = 90°), plotted in solid line. These parameters also agree with the timings 
of the primary and secondary eclipses, shown in dash lines. Clearly a pulse similar to that 
of KIC 6864859 can be detected without the aid of eclipses. Since the geometric bias can 
be well estimated, one could make consistency checks between the systems detected among 
the eclipsing binaries and the ones without the eclipses. The light curves of a few other 
candidates (KIC 4459068, 5790807, 4949194) with periods ranging from 25 d to 80 d are 
shown in Figure 4. KIC 4459068 and KIC 5790807 are associated with single eclipses (both 
of which are included in the EBcatalog by Slawson et al. (2011)) while KIC 4949194 does 
not exhibit any eclipses. Some other probable candidates include KIC 9172506 and 9028474, 
which show possible pulses with amplitudes similar to those of other variabilities in the light 
curve. Note that the eccentricity pulse signals have relatively large amplitude for early-type 
stars, raising the interesting possibility of studying the possible eccentric migrating binaries 
for those stars. 



5. Discussion 

The presence of more than a dozen highly-eccentric, long-period Kepler eclipsing binaries 
confirm that ~ 1% of stellar systems reside in the P-e plane near P(l — e 2 ) 3 / 2 = Pf ~ 15 d, 
as hinted by the Raghavan et al. (2010) sample. The binaries in this stream may be excited 
to high eccentricity in Kozai cycles and experience tidal migration in the process of forming 
close binaries. This interpretation has multiple implications that can be tested in the near 
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future: (1) A growing population at increasing periods is expected from the steady-state 
distribution, following dN/ dlgP oc P 1 / 3 . This can be probed by future releases of Kepler 
data, which will enable the exploration of longer-period regimes. One interesting possibility 
is the identification of the location from which these binaries originate. Larger samples 
can be obtained by using Radial- Velocity followup of EBs or by searching for eccentricity 
pulses in the Kepler data. (2) If affected by Kozai, these systems should have additional 
companions which can be searched for using high-resolution imaging. (3) Direct tests of 
Kozai oscillations can be achieved by preforming differential astrometry measurements of 
mutual inclinations between the inner binary and the companion (e.g., Muterspaugh et al. 
(2006)) for the nearby multiples in the stream. 
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log 10 P(d) 

Fig. 1. — Period and eccentricity of stellar systems presented in the Raghavan et al. (2010) 
volume-limited FGK dwarf multiplicity survey sample. Components within multiples are 
shown as filled circles while the binaries with no detected additional companions as crosses. 
The upper envelope in eccentricity for 10 day < P < 1000 day seems to be defined by a 
population of eccentric systems in multiples. They are located within a narrow range of 
constant angular momentum tracks at Pf = P{1 — e 2 ) 3 / 2 between 10 day and 30 day, plotted 
as solid lines. The period and minimum eccentricity of the long-period, highly-eccentric 
eclipsing binary sample described in this paper are shown in red arrows. The latter is 
robustly constrained by timing differences between primary and secondary eclipses tjj — tj. 
The dotted lines in red illustrate the observational cuts placed in searching for these eclipsing 
binaries. One is a lower cut in period of P > 40 d, and the other corresponds to an upper 
cut in eclipse timing difference tu — ti < 10 d, shown by the eccentricity lower limits that it 
sets as a function of period. Half of the duration for the current ( ~ 16— month), expected 
full (3.5 yr) and possibly extended (~ 10yr) Kepler data releases are plotted as arrows, 
respectively, at the top of the figure, illustrating roughly the expected completeness period 
threshold for eclipsing binary detections. 



- 12 - 



1 

0.95 

J as 

5 0.85 
0.8 
0.75 

1 

0.95 

I 

0.85 



E 



KIC 

P=li6d 

^L,., = 9.(M 



J 1 L 



J 1 1 L 



I 

! 

; 
i 
I 



_ 



J 1 1_ 



■ 

i ' i 



i 

-1 L 



200 



300 



400 



inn 



600 



0.8 



~7T 

;i:;:- : >--v '; 

5t LII -6.3d 



T 

: 
-: 

I 



_, , 1 — 



— i 1 1 1- 



_, , 1 1_ 



T 

■ t 

■: 



■'I 

■A 



I 
■ 

-. — 



1 - 



1 



0.98 



0.92 
0.9 



-i r 



300 



i — m 



500 



KIC 
■SEH71E 



— i 1 1 1 1 1 1 1 1 1 1— 



-i r 



600 



0.96 |Si liM =g.4d 



J 1 U 



■ 



— I 1 1 L- 



-i 1 1 u 



_1 1 L. 



_1 1 1_ 



200 



300 400 
BJD-2454833 



500 



600 



Fig. 2. — Light curves for the three longest-period eclipsing binaries in our sample. Primary 
and secondary eclipses can be readily recognized. The timing difference between the two 
eclipses place a robust lower limit to the orbital eccentricity, as discussed in section 4.1. 



13 



period=40.8778 



1 1 1 1 1 1 1 
1 1 


I 

:l 
•1 

# ! 
'•it 

V 

1 
1 
1 

1 1 


1 1 1 


i 

i 
i 
i 



12 14 16 18 20 22 




Fig. 3. — Folded light curves of KIC 6864859, a 41-d eclipsing binary that exhibits "eccen- 
tricity pulses" with ~ 2.5 x 10 -4 amplitudes between the two eclipses. Points with different 
colors represent observations in different period folds. The upper panel shows the portion 
of the light curve with two eclipses. The lower panel plots the same portion of the light 
curve with a zoomed view to show the eccentric pulses which are well described by the 
linear-perturbation model of ellipsoidal variations (Kumar et al. 1995) plotted in solid blue 
line with e = 0.63. The dashed lines indicate the eclipse timings required by the ellipsoidal 
variation model, which match well with the observed values. 
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Fig. 4. — Examples of systems exhibiting eccentricity pulses (see section 4.2). The upper 
panel shows a folded light curve (KIC 4459068, P=25 d) while the light curves shown in the 
middle and lower panels are not folded (KIC 5790807, P=80d and KIC 4949194, P=41d). 



